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ABSTRACT 

A set of coefficients has been obtained by 
fitting an empirical non-linear limb- darkening 
law of the fora 

“ 1 - Aj (l-'O-B^ • i log t 
to values of I(u)/I(l) computed from a grid of 
hydrogen line-blanketed model atmospheres. In 
addition, a sot of coefficients for the conventional 
linear law has been calculated. The comparison 
with previous theoretical values for tlxe linearized 
law shows good agreement with Grygar's (19<3o) results 
but systematically smaller vnlnes than those found 
by Ilosokawn (1967). 

! . INTR ODUCTION 

A number of detailed model atmospheres have been re-analyzed 
by Grygar (1965) to derive improved limb-darkening coefficients 
for early type stars. He noted that, although models from 
several different authors were used, good internal consistency 
existed among the separate determinations of the coefficients. 

A tabulation of mean coefficients is given for wavelengths 
ranging from ' 1000H to 900<)ft. This tabulation is restricted 
to surface gravities with log (g) =4.0 but additional 
coefficients from three of Strom's (1964) log(g) = 3.0 models 
are also depicted. More recently, Hosokawa (1967) has 
determined the darkening coefficients for a wider range of 
effective temperatures by analyzing the model atmospheres of 
Strom and Avrett, v 1966) for early type stars and of Glngerlch 


(19GG) lor cooler stars. Again the determinations arc only 
for log (g) = 4.0. This work differs from Grygar's in that 
tho Linearized darkening coefficient is defined by a flux 
equivalent relation rather than a least squares representation 
with unequal weights. w© shall adopt lfosokawa 's definition 
for and our results will l>e directly comparable. 

It is well known that except for the case of the simple 
gray atmosphere solution, the direct respresenta tion of tho 
emergent intensity distribution is markedly non-linear. To 
obtain observational confirmation of this relation by tho 
study of eclipsing binaries is difficult. Nevertlieleas , 

Grygar (1903) found that the photometric solution for the 
eclipsing system All Aurigae was improved, as indicated by 
a reduction in the residuals, when a non-linear darkening 
law was applied. This is not to say that solutions for all 
types of eclipsing binaries can bo improved by adopting a 
non-linear law. Most likely those systems which exhibit 
photometric (and spectrographic) complications, large outside 
eclipse variations, or are comprised of tldally or 
rotationally distorted components will not respond to this 
elaboration. Conversely, it cannot be expected that tliese 
systems will provide observational verification of the computed 
coefficients. 

In this present work a set of coefficients appropriate 
to the linear law of darkening (Hosokawa , 1967) and a sot 
for an eronirical non-linear law are calculated for early 
type stars. The latter set of coefficients is suitable for 
application in direct computer analysis of binary systems. 
Finally, the newly computed coefficients are compared both 
with previous tlieoretical values and with observed values. 

1 1 • MODEL CALCULATIONS 

A grid of flux corrected radiative and hydrostatic 
equilibrium model atmospheres with hydrogen line-blanketing 
has been used to compute I, (u) at wavelengths of astrophysical 
interest. The models cover tlie range 10000K * T ^ s 40000K; 
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2.5 log g 4.5 anti have a composition (by mnss) of X » 2/J, 

Y » l/J. flux ia const mt to within n few tenths of a 
percent. The opacity source** included aro the bound-free 
and the free-free transitions of hydrogen, helium and their 
ions both positive and negative (Mihalas, 1965; Vnrdya, 196.) , 

Gelt man, 1902; Fisc he 1, 1903 and Mclowell, et al. , 1966), 

electron scattering, and tho bound- bound transitions of the 
Lyman and Palmer series of hydrogen (Griera, 1960 and Underhill, 
1962). These models are currently being prepared for publication. 
No metal line blanketing has been included in these calculations 
lienee, the results in the far UV must be considered of lower 
weight than tliose in the near UV and visible. 

The emergent Intensity ns a function of (Chandrasekhar , 

1950) is defined as 

CD 

I ( , t ~ C)= S, (t )e‘ t A dt (1) 

0 ~7T 

where u is the direction cosine of the line of sight to the 
surface normal and S : (t, ) is tho monochromatic source function 
Including the scattering term (ICourganof f , 1963) . \ ten point 

Gauss i. a guerre quadrature formula was used to evaluate 
equation 1 in which x =* ty was taken as the independent 
variable. In order to verify the accuracy of tho calculations, 
values of i were chosen to permit a Gaussian integration of 
tlie moment equation ^ 

F (t - 0)- 2 1 ) < :) Um (2) 

which could then be compared with the net emergent flux, F, 
(Kourganoff, 1963) computed directly from 

CD 

( t n - 0) - 2 S (t ) K 0 (t )dt, (3) 

o " 

In all cases the difference between the emergent flux as 
computod by equations 2 and 3 is less than 1/3%. 

A linearised representation of tho emergent intensity 
as a function of Is commonly assumed. For this present 
work, llosokawa's (1967) relation defining the linear 
coefficient u^ is adopted. Tho coefficient is found from 
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tho equivalence relation 
1 K-.) 


icn 


a 


a - 
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which after Integration yields 
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In addition to the linearised loin, tho following empirical 
equations which express the intrinsic non- linearity of the 
emergent Intensity with u were assumed. The first can be 
recognized as that due to van't Veer (I960). 

Is I < )/I, (1) - 1 - A, (1- ) - B, <l-.i> 3 

II: I <")/I, (1) - l - A, (I- ) - B) u log 10 i 

'Flic coefficients, A, and li , v/ore found by the method of 
least squares. Probable errors ranged from 0.001 to 0.01 
with consistently smaller errors obtaining with the second 
empirical law. Tables 1 to 5 list tho values of Aj , 
pertaining to the socond law and the linear coefficients 
for each wavelength and each model. Separate evaluations 
were made at each surface gravity. Effective temperature 
varies across the table and wavelength varies down the 
table. 


III. D ISCU SS IOW 

It can be seen from the tabulated results that the limb 
darkening coefficients decrease smoothly with increasing 
surface gravity for all wavelengths. 

Tile magnitude of the effoct is rather small. On the other 
hand, the values vary markedly with temperoture and 
wave length in a manner similar to that found by previous 
investigators. It should be noted that the values of 
at ' =* 3362$ deviate from an apparent smooth wavelength relation. 
At this particular wavelength, tho wings of tho Halmer hydrogen 
lines alter the continuum opacity and produce tho nnamolous 
result. At tho present time, very few high quality limb • 
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darkening determinations for oarly type stars are available 
(e.g., refer to Wood (19632). The general tendency that the 
observed limb darkening coefficients are larger than the 
computed Is strengthened by those present results. 

Although the coefficients are not tabulated here, several 
models sere solved to obtain linearized Hmb darkening 
coefficients in the mannor described by Grygar. Excellent 
agreement between the two sets of coefficients was found 
indicating an equivalence for the respective model atmospheres 
used In the determination. However, a comparative study with 
Hosokawa's tabulation of mean limb darkening coefficients, 
where identical definitions for the linearisation lav were 
assumed, shows a systematic difference in the coefficients 
nearly Independent of wavelength. The sense of the difference 
is that the present coefficients are of the order of .04 units 
smaller than Hosokawa's. This difference is probably not 
significant for the Inverse observational problem. 
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